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Resumen 


El objetivo de este estudio consiste en validar mediante simulación CFD 
la efectividad de una nueva metodología de diseño conformada por un 
conjunto de correlaciones actualizadas, que permita el diseño de cada 
uno de los componentes de la MBT, con el fin de mejorar su eficiencia. 
En este estudio se realizó una investigación rigurosa de la literatura de 
la MBT, donde se determinaron los parámetros y las ecuaciones de 
diseño más influyentes en la eficiencia máxima. Finalmente, el diseño de 
la MBT se realiza con las correlaciones más relevantes encontradas en la 
literatura y el diseño de la MBT se valida mediante pruebas 
fluidodinámicas. Se concluye que la metodología propuesta para el 
diseño de la MBT puede alcanzar eficiencias hasta del 83 %, lo cual es 
satisfactorio para resolver la falta de métodos de diseño completos para 
el dimensionamiento de los diferentes componentes de la MBT (inyector, 
rotor y carcasa) de acuerdo con las condiciones de flujo del lugar de 
instalación. 


Palabras clave: turbomáquinas, pico-hidroeléctricas, eficiencia, CFD. 
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Abstract 


The aim of this study is to validate, by means of CFD simulation, the 
effectiveness of a new design methodology formed of a set of updated 
equations, which allows the design of each of the MBT components in 
order to improve its efficiency. In this study, a rigorous investigation of 
the MBT literature was carried out, where the most influential design 
parameters and equations in maximum efficiency were determined. 
Finally, the design of the MBT ¡is carried out with the most relevant 
equations found in the literature and the design of the MBT is validated 
by fluid-dynamic tests. It is concluded that the proposed methodology 
for the design of the MBT can reach efficiencies up to 83 %, which is 
satisfactory to solve the lack of complete design methods for the sizing 
of the different components of the MBT (nozzle, runner and housing), 
according to the flow conditions of the installation site. 
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Debido a la economía en el mercado competitivo, la modemización y el 
constante crecimiento de la población, varios estudios pronostican un 
gran aumento en el consumo de energía eléctrica para las próximas 
décadas, especialmente en los países desarrollados. Desde 2015, la 
energía hidroeléctrica constituye aproximadamente el 61 % del total de 
energía renovable mundial. De esto, las pequeñas y microcentrales 
hidroeléctricas aportan alrededor del 7 % (Tesfaye-Woldemariam, Lemu, 
8 Wang, 2018). Esta tecnología ofrece un gran potencial de explotación 
como altemativa de generación de energía a partir de fuentes no 
convencionales de energía renovable, particularmente en zonas no 
interconectadas (ZNI) de países en desarrollo, con un impacto reducido 
en el ecosistema, en comparación con las grandes centrales 
hidroeléctricas (Paish, 2002). Por otro lado, representan la tecnología de 
generación eléctrica menos costosa en relación con la energía solar y 
eólica (Organization, UNID, 2016). 


Para aprovechar el potencial hidroeléctrico disponible 
implementado en pequeñas centrales hidroeléctricas se requiere el uso 
de turbinas eficientes. Las tecnologías tradicionales de los sistemas de 
generación hidroeléctrica están integradas por turbinas de alta eficiencia 
como Francis, Pelton y Kaplan, y en menor proporción por MBT 
(Paredes-Gutiérrez, Palacio-Higuita, 8 Goméz-Gómez, 2008), a pesar de 
las ventajas comparativas en su simplicidad de diseño, bajo costo de 
fabricación y poca variación de eficiencia con cambios considerables en 
las condiciones de operación (Dragomirescu, 2016). 
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Contrario a los diferentes enfoques de optimización tanto 
numéricos como experimentales realizados en las últimas décadas 
(Olade, 1987), las metodologías de diseños de las MBT deben 
actualizarse (Tesfaye-Woldemariam et al., 2018) debido a que 
actualmente existen nuevos materiales con mejores propiedades físicas, 
equipos y técnicas de fabricación que permitan obtener componentes 
mejorados. Así, una buena metodología de diseño para dimensionar los 
componentes de la MBT tendrá un impacto significativo en el 
rendimiento, lo que aumentará la factibilidad de implementación como 
tecnología de generación para las condiciones adecuadas de sitio (Q y 
H). 


La MBT es una turbomáquina que aprovecha la energía cinética del 
fluido y la transforma en momento angular, el cual es convertido en 
energía eléctrica mediante un generador (Ceballos, Valencia, Zuluaga, 
Del-Rio, 8: García, 2017). Con base en las investigaciones reportadas en 
el estado del arte, se implementa un diseño de turbina similar 
compuesto por tres cuerpos: inyector, rotor y carcasa, como se muestra 
en la Figura 1(A). En la MBT, el fluido ingresa por el inyector, 
aumentando la velocidad del fluido y dirigiéndolo hacia la entrada del 
rotor, manteniendo el ángulo de ataque y la velocidad constantes a lo 
largo de la posición azimutal correspondiente a la interfaz entre la salida 
del inyector y la entrada del rotor. El flujo cruzado transfiere la energía 
en dos etapas del rotor, como se detalla en la Figura 1(B). En la primera 
y segunda etapa, la energía del fluido es entregada en un 70 y 30 %, 
respectivamente (Chiyembekezo S. Kaunda et a/., 2014; Ceballos etal., 
2017) 
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Primera etapa 


Segunda etapa 
B) 


Figura 1. MBT: A) vista isométrica; B) vista lateral del rotor. 


En la literatura se encuentran diferentes metodologías de diseño 
de cada uno de los componentes de la MBT, sobre todo para el inyector 
y el rotor, basadas en estudios tanto numéricos como experimentales. 
En cuanto a modificaciones del inyector, Adhikari y Wood (2017) 
realizaron un estudio del efecto de la geometría de la pared posterior del 
inyector en la eficiencia de la MBT mediante simulaciones 
computacionales, mostrando un aumento del 18 % en el rendimiento 
hidráulico de la turbina sobre el modelo sin la modificación. Luego, en 
un diseño posterior, Adhikari y Wood (2018) hicieron una investigación 
CFD en la que proponen un dispositivo electrónico que regula la entrada 
del inyector frente a cambios de flujo, con el fin de mantener constante 
la velocidad y el ángulo de ataque del fluido al entrar en el rotor, 
logrando eficiencias de hasta el 88 %. Finalmente, en cuanto al inyector, 


Rantererung, Tandiseno y Malissa (2019) diseñaron una MBT de 5 kW 
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para una ZNI; este diseño consiste en el uso de inyectores múltiples; 
busca determinar el número de inyectores necesarios para el mejor 
desempeño de la turbina. El modelo fue probado con 1, 2 y 3 
inyectores; se concluyó que el mejor comportamiento del MBT ocurre 


cuando hay tres inyectores generando 4.259 kW. 


Respecto a modificaciones en el rotor, Sammartano, Arico, 
Carravetta, Fecarotta y Tucciarelli (2013), así como  Chichkhede, 
Verma, Gaba y Bhowmick (2016) llevaron a cabo un análisis CFD 2D y 
3D sobre el diseño óptimo de la MBT con base en criterios de mejora 
fundamentados en parámetros geométricos que afectan de manera 
significativa la eficiencia de la turbina, como la posición, número, forma 
de los álabes y el ángulo de ataque del fluido; se obtuvo una eficiencia 
del 86 % utilizando 35 álabes y un ángulo de ataque de 220, Por otro 
lado, en la literatura se encontraron modificaciones en los álabes del 
rotor. Arellano-Vilchez (2015) efectuó un estudio utilizando técnicas 
CFD, examinando el comportamiento de los álabes del rotor con y sin 
bordes afilados, considerando también el grosor de los álabes; concluyó 
que el uso de álabes afilados mejora el rendimiento hidráulico de la 
MBT. Finalmente, con respecto a las modificaciones geométricas del 
rotor, la investigación realizada por Popescu (2017) analiza mediante 
CFD las bajas eficiencias de la MBT causadas por la interacción del fluido 
con el eje del rotor y la recirculación del flujo en régimen de baja 
rotación, determinando que, cuando el rotor trabaja sin eje, la eficiencia 


aumenta un 5 %. 
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Según la revisión del estado del arte existen diferentes ecuaciones 
para diseñar el inyector y el rotor de la MBT. Sin embargo, cada una de 
estas ecuaciones se realizó por separado, por lo que no es posible 
mostrar si al momento de concatenar estas ecuaciones para el diseño de 
una MBT se garantizarán los mismos resultados que se presentan por 
separado. Por lo tanto, es evidente la importancia de contar con una 
metodología de diseño actualizada para el dimensionamiento de cada 
uno de los elementos de la MBT, utilizando ecuaciones desarrolladas 
mediante metodologías experimentales y/o numéricas. El objetivo de 
este estudio es validar, mediante simulación CFD, la efectividad de una 
nueva metodología de diseño formada por un conjunto de ecuaciones 
actualizadas, que permita diseñar cada uno de los componentes del 
MBT, con el fin de mejorar su eficiencia. 


Metodología 


Ecuaciones de gobierno 
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La Figura 2 es una ilustración esquemática de la MBT, cuyo propósito es 


mostrar todos los parámetros geométricos que se consideraron en este 


estudio. 


Pared posterior 


ÓN inyector 


Primera etapa 


Flujo de entrada h 
——=> 0: 


Figura 2. Ilustración esquemática del diseño básico de una MBT. 
Fuente: adaptado de Adhikari y Wood (2017). 
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Para comenzar con el diseño de la MBT, las condiciones de 
operación deben tomarse como parámetros de entrada de acuerdo con 
las condiciones del sitio, Q y H; 20 l/s y 0.5 m, respectivamente. 
Posteriormente, se inicia el diseño del rotor, utilizando los diámetros 
externos (D.,,) del rotor de acuerdo con la relación caudal/altura de la 
Tabla 1 (Paz, Carrocci, Filho, £ Luna, 2007). Luego, el diámetro interno 
(Di,¿) se calcula mediante la Ecuación (1); la velocidad de flujo (v) en la 
entrada del rotor se determina mediante la Ecuación (2) (Adanta, 
Siswantara, 8 Prakoso, 2018). 


Tabla 1. Selección del diámetro del rotor. Fuente: Paz et a/. (2007). 


Q (m/s) Diámetro externo del rotor (mm) 
yJH(m) 


0.02236 - 0.04743 


0.07906 - 0.11068 


0.11068 - 0.15812 


A NAAA 


0.68 = Pine (1) 

veliaz-9g*H (2) 
Donde C, corresponde al factor de pérdida dentro del inyector, 

idealmente; C, se considera 1, mientras que experimentalmente 
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Sammartano, Morreale, Sinagra y Tucciarelli (2016) reportan valores 
iguales a 0.95. Posteriormente, para la selección del ángulo de ataque 
del fluido (w), diferentes estudios numéricos y experimentales proponen 
valores entre 12 y 22% (Chichkhede et al/., 2016; Wapnito, Siswantara, 
Adanta, 8: Prakoso, 2017). Luego, se determina el ángulo de velocidad 
relativa (f,) con respecto al rotor utilizando la Ecuación (3) dada por 
Sammartano et al. (2013): 


2tan(a)= tanta) (3) 


Basado en estudios realizados por Warjito et al. (2017) y 
Sammartano et al. (2013) se selecciona f,. Los ángulos de los álabes 


con respecto a la tangente de los diámetros interior y exterior deben ser 
Bi» Y Bop, respectivamente (Figura 3). 
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Diámetro exterior 


Figura 3. Geometría de los álabes de una MBT. Fuente: adaptado de 


Sammartano et al. (2013). 


Luego se calcula el radio de los álabes (r) y el ángulo de curvatura 


de los álabes (9) usando las ecuaciones (4) y (5), respectivamente 
(Olade, 1987): 


9) 
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$ = 2tan”* ls 5 E ,) (5) 
Tm lhsin(f,) 


Para calcular el ancho del rotor se debe obtener el número óptimo 


de palas (N,) mediante la Tabla 2. 


Tabla 2. Selección del número de álabes del rotor. Fuente: Paz et al. 
(2007). 


Diámetro externo del rotor (mm) Número de álabes 


Luego h, se calcula mediante la Ecuación (6), que corresponde a la 
altura del inyector a la entrada del rotor; donde 6, es la apertura del 
inyector a la entrada del rotor dada por Adhikari y Wood (2017), como 
se muestra en la Figura 4. Con estos datos es posible calcular el ancho 
del rotor mediante la Ecuación (7) (Olade, 1987): 


“o -=0.37 (6) 


Rext'0s 
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— TrDeypUsin(a)Za 


Donde z, es el número de álabes mojados en la primera etapa de 
la turbina, como se muestra en la Figura 2; el ancho del inyector se 


calcula mediante la Ecuación (8), y el régimen de rotación de la turbina 


(Uma, ) Se determina mediante la Ecuación (9) (Adhikari 8 Wood, 2017): 
_ Ww 
B= (8) 
Omax'Rext _ 1 nó 
v o 0 - 27) (2) 


s 


La forma geométrica de la pared posterior del inyector se puede 
obtener mediante la Ecuación (10) (Adhikari 8 Wood, 2017), donde 9 


tomará valores de O a 9, (Figura 4): 


Rio) =Ry + ho (1- $) (10) 


s 
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A 
y j Pared posterior 


ÓN inyecotr 


Flujo de entrada 
——> 


Figura 4. Ilustración esquemática de la pared posterior del inyector. 
Fuente: adaptado de Adhikari y Wood (2017). 


Finalmente, la eficiencia (n) de la MBT es determinada por la 
Ecuación (11) (Ceballos et a/., 2017), la cual corresponde a la relación 


entre la potencia de salida del eje (P,... =T :w) y la energía total del flujo 


out 


de agua (P;, =y*Q:H- 1): 


q= (11) 
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Modelo computacional 


Luego de realizar la metodología de diseño de la MBT con base en las 
ecuaciones de gobiemo encontradas en el estado del arte, se propone 
un diseño simplificado de los elementos (rotor, inyector y carcasa) que 
constituyen la turbina a partir de las ecuaciones (1) a (12); luego se 
realiza una operación booleana para obtener el volumen de control 
interno de la máquina hidráulica, el cual es mallado o discretizado en el 
módulo Mesh del software Ansys 2019R30, utilizando elementos 
tetraédricos con funciones avanzadas de curvatura y proximidad, y un 
tamaño de elemento mínimo de 2 mm, que fue suficiente para 
determinar una independencia de malla (Figura 5) con ocho millones de 


elementos con un error relativo menor al 1 %. 
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100 


Eficiencia (%) 


2 4 6 8 10 12 14 
Número de elementos (1e6) 


Figura 5. Independencia de malla. 


Con base en diferentes investigaciones donde se comparan los 
diversos modelos de turbulencia en problemas similares, y en algunos 
modelos computacionales distintos, se encuentra que el modelo 
numérico estándar k — e representa correctamente el flujo dentro de la 
MBT. Según diferentes estudios y guías de usuario de herramientas de 
simulación, este modelo es robusto, económico y proporciona una 
precisión razonable no sólo para el problema actual, sino también para 
una amplia gama de problemas. Por tanto, en este trabajo se utilizó el 
modelo estándar para análisis dependiente del tiempo, simulaciones 
fugaces con un tiempo de simulación de un segundo y con un paso de 


tiempo de 0.001 segundos, considerando el modelo bifásico homogéneo 
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compuesto por agua y aire a 25 *C, debido a que tiene un costo 
computacional menor en comparación al no homogéneo. La Figura 6(A) 
presenta la configuración CFX AnsysO; el fluido ingresa al sistema por el 
inyector (entrada) con una velocidad de 2.3 m/s; luego pasa por el 
dominio del anillo, donde se configuran las interfaces, con el fin de que 
pueda haber una transferencia de datos entre los dominios giratorio 
(rotor) y estacionario (carcasa e inyector). La primera interfaz está 
configurada entre el inyector, la superficie del anillo externo y la 
carcasa, como se muestra en la Figura 6(B). Para esto, se implementa 
una interfaz de tipo "fluido-fluido", con un método de conexión de malla 
definido por un modelo de interfaz de "conexión general" y sin marco de 
modelo de cambio/mezcla para dominios estáticos. La segunda interfaz 
está configurada entre la superficie del anillo interno y la superficie 
externa del rotor, como se muestra en la Figura 6(C). De manera 
similar, a la primera interfaz, se configura una interfaz "fluido-fluido" 
con "Estator de rotor transitorio" para el marco del modelo de 
cambio/mezcla debido a las condiciones de fluido transitorias entre el 
inyector y el rotor. 
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Pa 
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A) 


Figura 6. Configuración del CFX: A) condiciones de frontera; B) interfaz 
carcasa-inyector-anillo; C) interfaz rotor-anillo. 


Resultados 


En este apartado se presentan los resultados obtenidos mediante CFD, 


luego de implementar la configuración del modelo estudiado en CFX 
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Ansys6O) para una altura neta de 0.5 m y una simulación de flujo real de 
16.2 l/s debido a que la presión fue configurada como condición de 
frontera a la entrada. Se procesan y examinan los contomos de fracción 
volumétrica de agua, velocidad y presión obtenidos por el CFD, para 
caracterizar la influencia de la metodología propuesta en el desempeño 
de la turbina. 


En la Figura 7 se presenta la eficiencia con base en la relaciones 
de velocidad de la MBT para los resultados experimentales obtenidos por 
Sammartano, Morreale, Sinagra, Collura y Tucciarelli (2014), y los 
resultados numéricos obtenidos en este estudio, para el cual se 
realizaron seis simulaciones, variando el régimen de giro del rotor de a 
100 a 200 RPM, obteniendo una eficiencia máxima del 83 %, con un 
régimen de giro igual a 160 RPM, similar al obtenido mediante la 


Ecuación (9) propuesta por Adhikari y Wood (2017), con una relación de 


velocidad de (2) =17. Con base en tales resultados, este estudio 


muestra similitudes con los resultados experimentales obtenidos por 
Sammartano et al. (2014), quienes lograron eficiencias del 82.1% a 


relaciones de velocidad (2) cercanas a 1.7, como se muestra en la 


Figura 7. 
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O AKÁ 4 
penes V. Sammartano et al., 2014 
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Figura 7. Relación de velocidades de la MBT. 


En la Figura 8 se muestra la fracción volumétrica de agua en el 
plano de simetría (XY) de la MBT, donde se puede observar que existe 
una separación de flujo insignificante en la primera y segunda etapa 
debido a que el diseño implementado en el inyector permite una mayor 
similitud entre los ángulos de la velocidad relativa del fluido y la posición 
de los álabes (f, =f$,,) a lo largo de la posición del azimutal en la 
entrada del rotor; esto mejora las condiciones de flujo en la entrada y a 
través del rotor, lo que provoca una mejora en el rendimiento hidráulico, 
coincidiendo con los resultados obtenidos por Adhikari y Wood (2017). 


Además, se puede observar que el flujo no presenta alteraciones dentro 
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del rotor debido a que no hay un eje que lo atraviese. En consecuencia, 
el rendimiento hidráulico mejora en comparación con un rotor con eje, 
como concluye el estudio realizado por Sammartano et al. (2013). 


Fracción volumétrica de agua 


Plane 1 
0.999 


0.749 


0.500 


0.250 


0.000 


Figura 8. Fracción volumétrica de agua en la MBT. 


La Figura 9 presenta los vectores de velocidad trazados por la 
magnitud de la velocidad en el plano de simetría (XY). En este caso, se 


produce un aumento de la velocidad en la reducción de la sección 
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transversal del inyector, obteniendo una eficiencia del inyector 
determinada por la Ecuación (2) de alrededor del 97 %. Luego hay una 
reducción de la velocidad a la salida de la primera etapa debido a que el 
fluido entrega la mayor parte de su energía cinética en la primera etapa 
(Shepherd, 1956). Inmediatamente después se genera el flujo cruzado y 
el fluido transfiere el resto de su energía cinética a los álabes del rotor 


durante la segunda etapa. 


Vectores de velocidad 
3.0 


2.3 


0.8 


Figura 9. Vectores de velocidad del agua en la MBT. 


La Figura 10 presenta los contornos de presión del fluido en el 


mismo plano de simetría XY, donde se puede ver el paso del fluido a 
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través del MBT desde la entrada hasta la salida, donde los valores 
asignados a los colores verde y azul indican las presiones máxima y 
mínima generadas, respectivamente, por el chorro de agua. Es posible 
observar una pérdida de presión del fluido al pasar por la turbina, que 
se debe a dos factores: la conversión de H en energía cinética y las 
pérdidas en el inyector. Además, también se evidencian presiones de 
vacío que ocurren en los lados convexos de los álabes durante la 
primera y segunda etapa, siendo más evidentes en la primera etapa 
debido a que todos los álabes están sumergidos; esto concuerda con los 
perfiles de velocidad presentados en la Figura 9, donde se obtienen los 
mayores deltas de velocidad en los álabes ubicados en la primera etapa, 
lo cual es aconsejable debido a que la MBT funcionará como turbina de 
reacción, al utilizar la diferencia de presiones en los álabes para 
aumentar el momento angular del corredor (Ceballos et a/., 2017). 
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Contornos de presión 
6.0 


Figura 10. Contornos de presión manométrica en la MBT. 


Conclusiones 


Este artículo presenta la validación numérica de una metodología de 
diseño nueva y simple para el dimensionamiento de la MBT, que se 


utiliza por lo regular en todo el mundo en sistemas hidroeléctricos de 
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pequeña escala. Para tal efecto, se ha realizado un estudio analítico, 
utilizando diferentes ecuaciones para el diseño geométrico de cada 
elemento de la MBT, especialmente para el rotor y el inyector, a partir 
de investigaciones numéricas y experimentales que se utilizaron para 
establecer la metodología presentada, la cual fue validada mediante 
simulaciones tridimensionales de Navier-Stokes promediadas por el 
método de Reynolds (RANS) con el modelo de turbulencia k-e, un 
modelo de flujo de superficie libre homogéneo de dos fases, y el 
software comercial ANSYS CFX. 


La metodología establecida en este trabajo permite diseñar 
cada uno de los elementos de la MBT con base en las condiciones de 
operación del lugar de instalación (Q y H), logrando eficiencias 
hidráulicas del 83 % y acortando de forma significativa el tiempo 
requerido en la etapa de diseño. Se verifica que la metodología de 
diseño propuesta para el rotor y el inyector permite al diseñador obtener 
una MBT de alto rendimiento, siguiendo las condiciones de sitio (Q y H). 
Además, ofrece una alternativa actualizada a la metodología propuesta 
por Olade (1987) para el diseño de este tipo de turbinas. 


Trabajo futuro 
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En cuanto a trabajos futuros, conviene establecer diferentes ecuaciones 
para un nuevo diseño de álabes para el rotor en términos de las 
condiciones del sitio (H y Q), integrando nuevos materiales y técnicas de 
fabricación que permitan fabricar álabes con geometrías mejorada, que 
puedan incrementar la transferencia de energía en la primera y segunda 
etapa. Además, las modificaciones en la pared posterior del inyector 
podrían mejorar el perfil de velocidad en la entrada del rotor y la 
eficiencia de la turbina. 
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